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ABSTRACT 

NGC 2915 is a nearby blue compact dwarf with a differentially rotating H I disc ex- 
tending out to ~ 5 i?-band R25 radii. This disc serves as an ideal tracer of the system's 
gravitational potential in regions of the galaxy that are dominated by dark matter. 
We use new H I synthesis observations of NGC 2915, obtained with the Australia Tele- 
scope Compact Array, to search for non-circular flows within the outer Hi disc. Two 
independent methods are used, and the results of each interpreted in the context of 
relevant axisymmetric and non-axisymmetric perturbations of the potential. We find 
evidence for: (1) elliptical streaming associated with the spiral structure of the Hi disc 
and the central bar-like feature in the mass distribution, (2) a spherical dark matter 
halo, and (3) an axisymmetric radial outflow of ~ 5-17 km s~^ (~ 6-20 percent of the 
circular speed). A possible bar-like perturbation of the potential hinders attempts to 
unambiguously detect kinematic signatures of radial flows in the H I velocity field. The 
radial outflows are inconsistent with the plausible disc formation scenario in which gas 
from the surrounding inter-galactic medium is deposited on the outer H I disc and then 
transported towards the centre of the galaxy. They are, however, consistent with the 
possibility of some material being re-distributed towards the outer disc in order to 
conserve angular momentum as material flows inwards along a bar. 
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1 INTRODUCTION 

Observed velocity fields of late-type galaxies are often 
treated as being consistent with a purely circular flow pat- 
tern. However, the intrinsic, and therefore also the observed, 
properties of a galaxy's kinematics will be affected by the 
presence of non-circular velocity components within the or- 
bits of observable test particles. These non-circular motions 
are caused in various ways, and can be either systematic or 
random in nature. On small length scales (< 0.5 kpc) for ex- 
ample, stellar feedback can significantly disrupt gas dynam- 
ics, setting up expandi ng shells within a galaxy 's inter-stellar 
medium (ISM) (e.g. IWalter Brinksl IiQQqI I. On larger 
scales, stellar feedback can drive galac tic outflows in sys- 
tems of sufficiently low dynamical mass ( Meurer et al. ]|l992l: 
Martini [T996I : iMac Low k Ferraralll999l : Ivaii Evmeren et al.1 



20091 ). Globally, a non-axisymmetric gravitational potential 
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will force orbits to be intrinsically elliptical rather than cir- 
cular (Binnev 1978; Bosma 1978). 

Given an observed velocity field, it is the task of the 
investigator to determine what information the observed 
fiow pattern contains about the mass distribution and other 
properties of the system. In the case where the orbits of 
tracer particles are significantly affected by non-circular mo- 
tions, a careful characterisation of the non- circular signa- 
tures present within the veloci t y field is required. Several 
authors (e.g. iFranx et al. 19941: ISchoenmakers et all Il997l : 



ICanzian Allenlll997l : ISpekkens Sellwoodll2007l l have de 
veloped methods for investigating and interpreting complex, 
non-circular velocity fields of galaxies in the context of per- 
turbations of the potential. These methods allow further in- 
sight into the fundamental nature of the processes respon- 
sible for generating the non-circular components, as well 
as t he non-circular co mponents themselves. As an exam- 
ple, ^^^^^| |20Q4l) used the CO and Hi velocity fields 
of seven nearby spirals to search for evidence of radial gas 
fiows. They found deviations from pure circular rotation at 
the level of 20-60 km s~^ and attributed them to the ef- 
fects of bar streaming, infiows and warps. More recently, 
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Figure 1. H l surface density map of NGC 2915 from E2010a. The 
logarithmic intensity scale is in units of Mq pc~^. The solid el- 
lipse with semi-major axis length 98 arcsec represents the i?-band 
i?25 radius. The dashed ellipses have semi-major axis lengths of 
200 arcsec and 510 arcsec, and delimit the portion of the H I disc 
for which an harmonic decomposition of the line-of-sight velocities 
is carried out in this paper. 



iTrachternach et all (|2008l ) quantified the non-circular mo- 
tions in a sample of 19 gala xies from the Th e Hi Nearby 
Galaxy Survey (THINGS, Wa lter et al.ll2008l ). They found 
a median absolute amplitude of the non-circular motions, 
averaged over their entire sample, of 6.7 km s~^. From the 
measured non-circular velocity components they determined 
the mean elongation of the gravitational potential to be 
0.017 ± 0.020, consistent with a round potential. 

In this work we search for non- circular fl ows in the ex- 
tende d Hi disc of the nearby (4.1 ± 0.3 Mpc, iMeurer et al.l 
I2OO3I I galaxy, NGC 2915. Optically, the galaxy is classi- 
fied as a blue compact dwarf. Centre d on the o ptical disc 
{R25 ^ 98 arcsec in the i?-band, Meurer et al.l 1994) is a 
differentially rotating H i disc with well-defined spiral struc- 
ture, extending further than ^ 500 arcse c. We use new H i 
synthesis observations (jElson et al.ll20ld hereafter E2010a) 
to study the H I kinematics of the galaxy. The H I total inten- 
sity map of NGC 2915 is shown in FigfTJ In this paper we are 
primarily concerned with the search for and the physical in- 
terpretation of non-circular velocity components within the 
H I disc. Two independent methods are used to quantify the 
non-circular flows. The results are interpreted in the con- 
text of various axisymmetric and non-axisymmetric pertur- 
bations of th e potential. The first raethod uses the approach 
developed bv'Sp ekkens SellwoodI (|2007l l for fitting a gen- 
eral axisymmetric or non-axisymmetric model to the veloc- 
ity field of a galaxy. They develop the model specifically 
for the case of bi-symmetric distortions (bar- or oval-like) 
and for purely axisymmetri c radial flo ws. The second tech- 
nique is that of Franx et al. (1994) and Schoenmakers et al. 
(|l997l ). used to interpret the harmonic decomposition results 
of an observed velocity field in the context of time- varying 
perturbations to an axisymmetric potential. 

The outline of this paper is as follows: the H i data are 
introduced in Section [2l The observed H i velocity field of 
NGC 2915, along with a pure circular rotation model ap- 
proximation of it, is shown in Section [S] We proceed to 
decompose and model the observed line-of-sight velocities 
in Section O using each of the two above-mentioned inde- 
pendent methods. The various results are then interpreted. 



compared and discussed. Section [4.2.41 focuses on the possi- 
ble tri- axial nature of NGC 2915's dark matter halo. Finally, 
our summary and conclusions are presented in Section [S] 



2 H I DATA 

In this work we utilise the Hl-line data set of NGC 2915 
presented in E2010a. The Hi observations were carried out 
between 2006 October 23 and 2007 June 2 (project number 
C1629) with six different Australia Telescope Compact Ar- 
ray antenna configurations using all six antennas. E2010a 
describe the various telescope setups of the observations as 
well as the reduction procedures used to arrive at the cali- 
brated, continuum-subtracted, deconvolved Hi data cubes. 
This work is based on their naturally- weighted Hi data 
cube which has spatial and spectral resolutions of 17.0 arc- 
sec X 18.2 arcsec and 3.2 km s~^ respect iveljO. We use the 
Hi velocity field which was extracted from the naturally- 
weighted Hi data cube by fitting a third-order Gauss- 
Hermite (/is-GH) p olynomial to each Hi line profile. The 
reader i s referred to Franx et al.l (| 19941 ) ; iNoordermeer et aD 
(|2007[ ): Ide Blok et al.1 (|2008l l for various discussions of the 
merits of such a means of velocity field construction. In the 
context of this work it suffices to say that the /is-GH param- 
eterisation of a line profile is less susceptible to the effects of 
profile asymmetries on circular velocity estimates made for 
a rotating disk. 



3 EVIDENCE OF NON-CIRCULAR FLOWS 

The /i3-GH H i velocity field of NGC 2915 is shown in the top 
panel of Fig. [2] The general pattern of the iso- velocity con- 
tours is that of a disc that is rotating linearly at inner radii 
and differentially at outer radii. The sharp kinks and wiggles 
along the iso- velocity contours, particularly prevalent at in- 
ner radii, are indicative of non-circular velocity components 
within the Hi. E2010a fitted two different tilted ring mod- 
els to this H I velocity field. The only significant difference 
between their two models was the best-fitting inclination 
profile. Their so-called model CI had a constant inclination 
of ^ 55° for all rings, while their model SI exhibited a sharp 
increase of ^ 20° when moving from the outer to the in- 
ner disc. The models differed only slightly in terms of their 
circular velocity profiles (see E2010a for the full details and 
results of the modeling procedure). 

The middle panel of Fig.[2]shows a purely circular model 
velocity field generated using the best-fitting parameters of 
tilted ring model CI from E2010a. This model cannot ade- 
quately account for the observed velocities. As a measure of 
the typical departures from pure circular motion at various 
positions within the galaxy, a residual velocity field (shown 
in the bottom panel of Fig. [2]) is produced by determin- 
ing the absolute difference between the observed and model 
velocity fields. From this residual map it is clear that the 
vast majority of the large residuals lie close to the centre of 
NGC 2915 and along the Hi spiral arms. The mean absolute 



At the distance of NGC 2915 (~ 4.1 Mpc), 1 arcsec corresponds 
to a linear distance of 19.8 pc. 
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Figure 2. Top panel: /is-GH velocity field from E2010a; Middle 
panel: model circular velocity field based on the tilted ring model 
from E2010a. Bottom panel: absolute residual velocity field ob- 
tained by subtracting the model field from the /is-GH velocity 
field. In all panels, contours are separated by 15 km s"-*^, with the 
thick contour marking the systemic velocity at 465 km s"-*^. Grey 
and black contours represent receding and approaching halves of 
the galaxy respectively. The intensity scale of the line-of-sight ve- 
locities is in units of km s"-*^. All panels show the same area. The 
white ellipse represents the i?-band i?25 radius of 98 arcsec of the 
stellar disc. The white cross represents the photometric centre 
determined my E2010a. 



residual is '-^ 7.8 km s~^, with 69 percent of the area of the 
residual velocity field having a value less than 10 km s~^ 
(the approximate velocity dispersion of the outer H i disc) . 
This suggests that a significant fraction of the Hi is regu- 
larly rotating in circular orbits about the kinematic centre. 
However, almost 9 percent of the residuals are larger than 
15 km s~^, some being as high as 20-25 km s~^, indicative 
of significant non-circular H I velocity components. 



4 QUANTIFYING THE NON-CIRCULAR 
FLOWS 

4.1 Method 1: velfit 

4.1.1 Introduction 

The first method used to quantify the non-circular flows 
within the H i disc of NGC 2915 is t he pr ocedure initially 
introduced by ISpekkens SellwoodI (|2QQ7f ). Given certain 
assumptions regarding the nature of the non-circular flows, 
their fitting routine, velfit, is able to fit to an observed 
velocity field a kinematic model describing the non-circular 
velocities. The line-of-sight velocity, V^os, of a test particle 
orbiting the kinematic centre of a galaxy within the plane 
of the disc is 

Vios = Vsys + sini(yt cos^ + Vrsin^), (1) 

where Vsys is the systemic velocity of the galaxy, Vt and Vr 
are the tangential and radial velocity components respec- 
tively, and is the angle in the disc plane relative to the 
projected major axis. As explained bv lSpekkens SellwoodI 
the above expression can be re- written in terms of a 
Fourier series: 

00 

Vios = Vsys + sin i[Vt cos + ^ Vm,t COS cos{mO + Om,t) 

m — 1 

00 

+ Vr sin 6 + Vm,r slu 6 cos{m6 + 6m,r)] , (2) 

m — 1 

where Vm,t and Vm,r represent the magnitudes of the non- 
circular velocity components, and Om,t and Om,r their phases 
relative to some axis. 

VELFIT makes the following assumptions in order to 
model Vm,t and Vm,r for a two-dimensional velocity field: 

(i) The non-circular motions are caused by a bar- or 
oval-like distortion to an otherwise axisymmetric potential. 
Alternatively, the velocity field is treated as containing a 
radial flow component. The harmonic order of the poten- 
tial perturbation in the plane of the disc is represented by m. 

(ii) The m — 2 terms dominate the non-circular velocity 
components for a strong bi-symmetric potential. VELFIT 
therefore only determines the V2,t and V2,r terms in 
Eqn. [2] VELFIT assumes that the bi-symmetric distortion 
has a fixed position angle, thereby yielding the fitted 
model insensitive to spiral structure. Alternatively, velfit 
can ignore all harmonic terms with m > 0, thereby fitting 
a purely axisymmetric radial flow model to the velocity field. 

(iii) The disc under consideration is not warped. 

Assuming the m — 2 term to dominate the non-circular ve- 
locity components, VELFIT fits the following model velocity 
field to the observed velocity field: 

Vmodel = Vsys + siu i[Vt COS - V2,t COs(2^fe) COS 6 

-y2,rsin(26>6)sin6'], (3) 

where Oh is the major axis of the bi-symmetric distortion 
(e.g. a bar) relative to the projected major axis. The rou- 
tine determines radial profiles for Vt, V2,t and V2,r that allow 
for the best fit of the model vel ocity field to the obser v ed ve - 
locity field (see Section 3.4 of ISpekkens k SellwoodI (|200/t ) 
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for a detailed discussion of the method by which velfit 
determines the radial profiles). The above-mentioned bi- 
symmetric model is also described by various disc param- 
eters, the values of which velfit estimates when fitting the 
model to the observed velocity field. These additional pa- 
rameters are: (1) Xc and ^c, the coordinates of the kinematic 
centre of the disc; (2) Vsys] (3) e^, the ellipticity of the disc; 
(4) (j)'^^ the major axis position angle of the disc in the sky 
plane; and (5) the bar position angle in the sky plane. 
Assuming all harmonic terms with m > to be negligible, 
VELFIT can also fit a purely axisymmetric radial flow model 
to the velocity field: 

Vmodei = Vsys + sini(Vt cos^ + Vr sin^). (4) 

Spekkens SellwoodI (| 20071) , as well as 

ISellwood Zanmar Sanched ('2OIO'), highlight several 
features of the velfit routine which they argue makes 
it better than other methods of estimating non-circular 
motions in rotating discs. VELFIT fits the entire velocity 
field at once, thereby allowing the effects of a mild bi- 
symmetric distortion to the potential to be detected in 
the flow pattern. Weak kinematic signatures, they say, can 
be missed when fitting models to only a small portion of 
the velocity field. They suggest further that considering 
only a bi-symmetric model allows for a less ambiguous 
identification of m = 2 components within the observed 
velocity field. Most importantly, they highlight the fact 
that their fitting routine does not rely on the assumption 
that the non-circular velocity components are small relative 
to the mean circular motion. 



17 arcsec and Rmax = 300 arcsec, softening parameters of 
Al = A2 = 1.5 X 10~^ were used to slightly smooth the radial 
variations of the Vt , Vr , V2 t and V2,r profiles of the hi^h- 
resolution H i velocity fielcfn. ISellwood Zanmar Sanched 
(2010 ) point out that the penalising effects of the softening 
parameters are insensitive to variations of a factor of a few 
about their values. 

Models were fitted to the entire low-resolution Hi ve- 
locity field, with no smoothing parameters applied. For all 
models, the differences between the observed and modelled 
line-of-sight velocities were weighted by the inverse of the 
typical velocity dispersion of the outer Hi disc (10 km s~^; 
E2010a). The position angle in the sky plane of the bi- 
symmetric distortion (i.e. the Hi bar, clearly visible in the 
Hi total intensity map) was not allowed to freely vary in 
the bi-symmetric flow models, and was fixed to = 322°. 
This value corresponds to the dynamical estimate of the po- 
sition angle for the central H i disc produced by E2010a who 
fitted tilted ring models to the high-resolution Hi velocity 
field. Our best-fitting bi-symmetric flow models were found 
to be virtually identical for variations of a few degrees of ^5. 
Non-parametric estimates of the true uncertainties in the es- 
timated parameters were generated by applying a bootstrap 
method. The method uses the spread of estimated parame- 
ter values from repeated fits to the observed H i velocity field 
to yield the non-parametric estimates. A full description of 
this method of error estimati on in the context of velfit 
i s pro vided in Appendix B of ISellwood k, Zanmar Sanched 
(|2010h . 



4-1.2 Fitted models 

VELFIT was used to fit bi-symmetric and radial flow model 
velocity fields (Eqns. [3] and [H respectively) to the /is-GH 
Hi velocity field of NGC 2915. Models were applied to 
the 17 arcsec-resolution velocity field (Fig. [3|), as well as 
a 34 arcsec-resolution velocity field (Fig. For all models, 
the parameters Xc, l/c, Vsys^ Cd, 0d and 0^ were allowed to 
vary freely. VELFIT allows the user to specify a galactocen- 
tric radius within which no pixels of the velocity field are 
fitted. In this way a central region of the velocity field can 
be excluded from being modelled. We, however, always ap- 
plied our models from the very centres of the velocity fields. 
Fitted pixels were sampled in ring widths of 17 arcsec. 

Stable fits to the high-resolution H I velocity field could 
not be obtained when modeling the entire disc. The resulting 
Vt, Vr, V2,t and V2,r profiles exhibited large, non-systematic 
variations at radii R > 300 arcsec. The high-resolution ve- 
locity field was therefore only fitted out to R — 300 arcsec. 
ISellwood Zanmar SanchezI (|2010l ) describe in their Ap- 
pendix A an optional penalty that can be applied within 
VELFIT to slightly smooth the radial variations of the fit- 
ted Vfc, Vr, V2,t and V2,r profiles, thereby allowing for more 
stable modeling results. They state 

0002 x^1, 

\ AR ) 

as an appropriate softening parameter, where Vtyp is an es- 
timate of the typical orbital velocity in the disc, AR is the 
ring spacing and Rmax is the maximum galactocentric ra- 
dius of the fitted model. Using Vtyp = 85 km s~^, AR = 



4.1.3 Results 

The final bi-symmetric and radial fiow models of the high- 
and low-resolution observed H i velocity fields are shown in 
Figs. [3] and |4] respectively. Also included in the figures are 
the corresponding absolute residual velocity fields, obtained 
by subtracting the models from the observed fields. The 
best-fitting Xc, yc, Vsys, (t>'d and parameters for each 
of the final models are presented in Table [1] The Xc and 
He values represent the respective absolute x and y offsets 
of the fitted kinematic centres from the photometric centre 
of ^2000 = 09^ 26"^ 12.6^ ^2000 = -76° 37' 37.8'', deter- 
mined by E2010a. The Vt, V2,t, V2,r and Vr radial profiles 
for the bi-symmetric and radial fiow models are shown for 
the high- and low- resolution H i velocity fields in Fig. (5] For 
comparison and for reference between plots, circular velocity 
profiles from tilted ring models fitted to each velocity field 
are presented as solid black curves in the various panels. The 
circular velocity profile used for the high-resolution velocity 
field is from the so-called tilted ring model CI presented in 
E2010a. The best-fitting geometric and orientation parame- 
ters of this model were assumed for a tilted-ring model fitted 
to the low-resolution Hi velocity field, the circular velocity 
profile for which is shown in the lower panels of Fig. [5] 

The overall kinematics of the models are consistent with 
those of the observations. The iso- velocity contours of the bi- 
symmetric fiow models of both the high- and low-resolution 

"^^ The magnitude of the penalty can be set independently for 
the axisymmetric and non-axisymmetric terms, we used the same 
penalty for both. 
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Figure 3. Kinematic models of the high-resolution NGC 2915 Hi velocity field. Panels a, b, and c show the observed velocity field, the 
best-fitting axisymmetric radial flow model and the best-fltting bi-symmetric flow model respectively. The ellipses in the top left panel 
delimit each of the model flelds. The velocity flelds are all plotted using the same intensity scale in units of km s"-*^. Iso- velocity contours 
are spaced by 15 km s"-*^, with the thick contour corresponding to 465 km s"-*^. Panels d and e show the respective residual maps. Both 
maps are plotted using the same intensity scale in units of km s"-*^. Kinematic models were only fltted to radii R < 300 arcsec, all maps 
therefore represent only this inner portion of the galaxy. The portion of the observed velocity fleld shown in the top left panel is delimited 
by a black rectangle in Fig. [4^. 




Figure 4. Kinematic models of the low-resolution NGC 2915 Hi velocity fleld. Panels a, b, and c show the observed velocity fleld, 
the best-fltting axisymmetric radial flow model and the best-fltting bi-symmetric flow model respectively. The ellipses in the top left 
panel delimit each of the model flelds. The velocity flelds are all plotted using the same intensity scale in units of km s"-*^. Iso- velocity 
contours are spaced by 15 km s"-*^, with the thick contour corresponding to 465 km s~^. Panels d and e show the respective residual 
maps. Both maps are plotted using the same intensity scale in units of km s"-*^. Note the difference in scale between these maps and the 
maps presented in Fig. [3] above. The black rectangle shown in panel a delimits the portion of the high-resolution Hi velocity fleld shown 
Fig. a. 



H I velocity fields are sharper and more twisted at innermost 
radii than those of the corresponding radial flow models. De- 
spite being kinematically simpler than the bi-symmetric flow 
models, it is the radial flow models that best match the ob- 
servations. This is particularly evident in the case of the low- 
resolution H I velocity field for which the radial flow model 
neatly reproduces the large-scale kinematic features. The 



residual patterns for the high-resolution Hi velocity fields 
are similar for both models. Residuals of ^ 20 km s~^ are 
observed at inner radii, and residuals < 10 km s~^ through- 
out the rest of the disc. These smaller residuals are clumped 
together on roughly 80-100 arcsec angular scales 1.6- 
2.0 kpc). Quite noticeably, the radial flow model of the low- 
resolution H I velocity field allows for small residuals at inner 
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Table 1. Summary of best-fitting model parameters for NGC 2915." 



1 2 3 4 5 6789 10 

Model Xc yc Vsys i (f^^ (f^b ^red ^ < l^^l > 

(arcsec) (arcsec) (km s"-*-) (deg) (deg) (deg) (km s"-*-) 



High res. 

Bi-symmetric 1.1 ±3.6 4.5 ± 1.5 462.4 ± 1.2 53.6 ± 1.1 304.4 ± 2.2 322.0 1.80 514 6.5 

Radial 1.9 ±7.1 4.1 ±3.7 462.2 ±2.1 53.5 ±3.1 295.6 ± 3.3 ... 1.94 533 6.1 



Low res. 

Bi-symmetric -5.2 ± 10.3 -20.5 ± 5.5 458.5 ± 1.7 49.1 ±3.0 293.3 ± 3.5 322.0 0.71 1177 8.5 
Radial 3.7 ±5.5 2.0 ± 4.3 460.3 ± 0.8 52.7 ±2.0 292.6 ± 2.2 ... 0.67 1202 3.7 



"Column 1: type of model velocity field; Column 2/3: absolute x/y offset of kinematic centre from the photometric centre; Column 4: 
systemic velocity; Column 5: disc inclination; Column 6/7: position angle of disc/bi-symmetric distortion; Column 8: minimum reduced 

statistic; Column 9: number of degrees of freedom when fitting; Column 10: average absolute residual. 
Model parameters best-fitted to the high/low-resolution Hi velocity fields are listed in the top/bottom sections of the table. Note that 
0^ is held fixed when fitting the bi-symmetric models, the only parameter with this distinction. 




J I U < \ < \ I I < \ < L 



100 200 300 400 500 100 200 300 400 500 
R (arcsec) R (arcsec) 

Figure 5. Fitted velocity components for NGC 2915 Hi velocity field. In each panel the solid black curve represents the circular velocity 
profile of a tilted ring model fitted to the Hi velocity field. The left-most panels show the Vt, V2,t, ^2,r components (red circles, green 
squares, and blue triangles respectively) of the best-fitting bi-symmetric flow models for the high- and low-resolution H I velocity fields 
(top and bottom panels, filled and open symbols respectively). The right-most panels show the Vt and Vr components (red circles and 
stars respectively) of the best-fitting axisymmetric radial flow models for the high- and low-resolution H I velocity fields (top and bottom 
panels, filled and open symbols respectively). In all panels, the error bars represent non-parametric estimates of the true uncertainties 
in the parameters, generated by applying a bootstrap method. 



radii, as well as small residuals that are more uniformly dis- 
tributed throughout the entire disc. This is partly due the 
the lower angular resolution of the data, yet also indicates 
the accuracy of the fitted model. Indeed, the average abso- 
lute residual corresponding to this model is only 3.7 km s~^. 
This model is also able to reproduce some of the closed iso- 



velocity contours seen in the north-west portion of the H I 
disc. 

The best-fitting disc parameters for both the high- and 
low-resolution Hi velocity fields are very similar to one 
another. Furthermore, the fitted disc parameters are sim- 
ilar to those obtained by E2010A who fitted a tilted ring 
model to the high-resolution Hi velocity field. The best- 
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fitting inclinations of 53° closely match their almost con- 
stant inclination of 55°, while the systemic velocity es- 
timates of 460 — 462 km s~^ agree well with their pre- 
ferred 465 km s~^. The fitted kinematic centres for the 
high-resolution H i velocity fields are within ^ 4.6 arcsec 
(~ 91 pc) of the photometric centre. The best-fitting posi- 
tion angles determined by E2010A varied from ^ 315° to 
^ 285° in a roughly linear fashion when moving from the 
inner to the outer Hi disc. The best-fitting position angles 
of all fitted models lie within this range. 

The radial variations of the mean orbital speed of the 
gas as derived using velfit (Vt, filled and open red circles 
in Fig. [5]) agree well with the corresponding circular veloc- 
ity profiles from E2010A (solid black curves in Fig. [5])- All 
of the final velfit models contain significant non-circular 
velocity components. The magnitudes of these components 
are similar for both the bi-symmetric and axisymmetric 
models, exhibiting similar slight radial variations. The ax- 
isymmetric models suggest the presence of radial flows of 
10 — 15 km s~^ for the region < R % 100 arcsec. 
Beyond i? ^ 100 arcsec, the inferred radial velocities in- 
crease in magnitude, reaching typical values of ^ 15 km s~^ 
for 120 arcsec < R < 250 arcsec. The outermost disc 
(R > 250 arcsec) contains radial non-circular motions with 
magnitudes of a few km s~^. The bi-symmetric models fit- 
ted to both the high- and low-resolution Hi velocity fields 
suggest the presence of bi-symmetric flows about an axis 
inclined ^ 18° — 30° relative to the kinematic major axis. 
The magnitudes of the V2,t and V2,r components are similar 
to one another at all radii, increasing from ^ 8 km s~^ 
at < 100 arcsec to - 15 - 20 km s"^ for 100 arc- 
sec < R < 230 arcsec. 

4- 1-4 Discussion 

For discussion purposes, we consider only the smoothly- 
varying radial profiles of the non-circular velocity compo- 
nents fitted to the low- resolution Hi velocity field. These 
profiles closely match the corresponding profiles of the high- 
resolution H I velocity field. Regarding the modeling results, 
an important question to answer is the following: given that 
the V2,t V2,r and Vr profiles vary with galactocentric radius 
in an approximately similar manner, should the optimal ra- 
dial flow model be preferred over the optimal bi-symmetric 

model, or vice versa? 

Using VELFIT, [Spekkens SellwoodI (|2007l l found Vr of 
the order of 10 km s~^ throughout the inner < 1.5 kpc 
(- 85 arcsec) portion of the Hi disc of NGC 2976. They 
argue that such large radial flows will lead to all of the de- 
tected gas in the quiescent system being displaced on kpc 
length-scales on time-scales of 1-3 Gyr. They therefore dis- 
favour the results of their optimal radial flow model. While 
it is true that significant radial flows within a gaseous disc 
will redistribute the gas over extended time periods, it can 
be argued that such radial flows are expected near the star- 
forming cores of galaxies. Indeed, several authors have linked 
the effects of stellar feedback to ra dially expandin g gas com- 
ponents in nearby galaxies (e.g. Ivan Ewneren et al. 200i; 
lYoung et al]l2003l : IWalter Brinkslll999l ). For the case of 
NGC 2915, we therefore do not exclude the possibility of ra- 
dial flows within the H i disc, especially at inner radii where 
the effects of star formation may dominate the gas kinemat- 



ics. In Elson et al., 2010 (in prep.) we produce detailed three- 
dimensional models of the Hi observations of NGC 2915 in 
order to link the central H i kinematics to the observed star 
formation activity. 

The H I total intensity map of NGC 2915 (Fig.p clearly 
displays a bar-like central Hi feature, thereby providing 
strong evidence for an m = 2 perturbation to the poten- 
tial of the galaxy. The stellar disc of NGC 2915 coincides 
in position, size and orientation to the observed bar-like H i 
feature. It can therefore be argued that the V2,t and V2,r ra- 
dial profiles of the bi-symmetric flow models are reasonable 
characterisations of Hi streaming motions of 0-20 km s~^, 
driven by a bar in the mass distribution. That the V2,t and 
V2,r amplitudes are relatively small for R < 100 arcsec is not 
surprising given the few resolution elements across the bar in 
this region. A slightly more subtle, yet arguably more perti- 
nent question to answer, is whether having fixed the position 
angle of the bi-symmetric distortion to = 322°, i.e. to that 
of the bar-like feature seen in the H i distribution, leaves the 
bi-symmetric models insensitive to bi-symmetric distortions 
that are not aligned with this feature. To address this con- 
cern, additional bi-symmetric flow models were fitted to the 
low- resolution H i velocity field, this time with allowed to 
freely vary. The models were fitted separately to the entire 
Hi disc as well as to only the outer Hi disc {R > 200 arc- 
sec). Modelling the entire disc yielded V2,t and V2,r radial 
profiles very similar to those presented in Fig. [5]:, with the 
corresponding V2,t and V2,r profiles differing respectively by 
3.4 km s~^ and 3.6 km s~^ on average for the cases in which 
05 was fixed and allowed to freely vary. The best-fitting po- 
sition angle in the plane of the sky for the bi-symmetric 
distortion is = 319.2°. The model fitted only to the outer 
disc has a best-fitting position angle of 0^ = 311.9°. These 
fitted position angles, being very similar to the position an- 
gle 05 = 322° of the observed bar- like feature in the mass 
distribution, suggest that this feature is indeed driving the 
elliptical streaming motions of the gas, and that our fitted 
models are not insensitive to bi-symmetric distortions that 
are misaligned wit h it. 

Spekkens & Sel lwoodI (|2007h . using velfit, found evi- 
dence for bi-symmetric flows of similar magni tude (^ ~ 
20 km s"^) for the case of NGC 2976. Masset Bureaul 
(2003) used hydrodynamical simulations to show that 
NGC 2915's observed Hi spiral structure can be accounted 
for by either an unseen bar or a rotating tri-axial dark mat- 
ter halo. A problem with the bar scenario is, however, that 
the mass of the required bar (~ 5 x 10^ M©) is very large rel- 
ative to the total stellar mass of the galaxy (- 8.0x10^ M©), 
thereby making the nature of such a bar problematic. 

4.2 Method 2: reswri 

Our second independent attempt at quantifying the non- 
circular gas motions in the Hi disc of NGC 2915 involves 
an harmonic decomposition of the line-of-sight-velocities 
and then, u s ing t he the oretical framework l a id do wn by 
iFranx et aD (|l994l l and |Schoenmakers et al." (1997), con- 
necting the measured kinematics to perturbations of the 
potential. These authors assumed a small stationary dis- 
tortion of the potential and then used epicycle theory to 
analytically determine the equations of motion of a test par- 
ticle. The equations of motion were used to determine the 
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result ing li ne-of-si^ht velocities. The m ethod of iFranx et alJ 
(|l994l l and ISchoenmakers et al ] (|l997l ) is therefore expected 
to yield reliable results only for small departures from circu- 
lar motion. In this sense, RESWRI is not as widely applicable 
as the VELFIT method which can fit arbitrarily large non- 
circular velocity components. It is, however, able to accom- 
modate warped discs with non-constant major axis position 
angles. RESWRI does not assume a particular type of poten- 
tial distortion, unlike the velfit technique. 



4-2.1 Fitted models 

Following ISchoenmakers et al the line-of-sight ve- 

locities are expressed as 



Vlos — Vsi 



' COS m -\- Sn 



(6) 



where N is the order of the fit, c^/ and s^/ represent the 
magnitudes of the non-circular line-of-sight velocity compo- 
nents, and is the azimuthal angle in the plane of the galaxy. 
Very importantly, because Eqn.[6]is a Fourier series for the 
line-of-sight velocities, m represents the harmonic order of 
the potential perturbation in the plane of the sky. This is 
unlike m in Eqn.[2]which represents the potential perturba- 
tion in the plane of the disc. The two harmo nic numbers ar e 
related by — mdil. Following iTrachternach et al.l (|2008l l. 
the line-of-sight velocities were expanded up to order = 3. 

The GIPSY task RESWRI was used to fit Eqn. [6]to con- 
centric rings in the high-resolution Hi velocity field. Unlike 
VELFIT, RESWRI does not use the entire observed Hi velocity 
field in a single fit, but rather divides the disc into concen- 
tric rings and then separately determines the non-circular 
line-of-sight velocities for each of them. As part of the de- 
composition process, RESWRI first fits a circular tilted ring 
model to the velocity field to obtain a measure of the ori- 
entation parameters of the rings as well as the kinematic 
centre for each ring. RESWRI then subtracts the circular 
rotation model from the observed velocity field and deter- 
mines the non- circ ular c omponents of the residual fi eld. Both 
iFranx et al.l (| 19941 ) and ISchoenmakers et al caution 
that the centre position of the tilted ring model should al- 
ways be kept fixed when running RESWRI. If allowed to vary, 
the centre will drift so as to make real C2 and S2 harmonics 
disappear. We therefore fixed the centres of the rings to the 
position of the photometric centre mentioned in Sec. 14.1.31 
The tilted ring model orientation parameters were allowed 
to vary freely, since incorrectly fixing them would result in 
incorrect measures of the harmonic components. The default 
uniform weighting was used during the fitting procedure to 
ensure that points along the minor axis contributed equally 
to the measured harmonics. These points contain the most 
information about non-circular motions. Ring widths equal 
to the full width at half maximum of the synthesised beam 
were used to ensure that adjacent rings were independent of 
one another. 

Rings with radius R < 200 arcsec were not included 
in the harmonic decomposition. The complex gas dynamics 
within this region are most likely dictated by the ongoing 
high-mass star formation at the centre of NGC 2915, and 
therefore cannot be meaningfully interpreted in the context 
of perturbations of an axisymmetric potential. Furthermore, 



the magnitudes of these non-circular motions are expected 
to be significan t fractions of the cor r espon ding circular ve- 
locities. Finallv. Schoenmakers et al.l (|l997 h caution that er- 
rors in the tilted ring fits are expected to be large in regions 
where V{R) oc thereby preventing a reliable measure of 
the harmonic components of the velocity field. For these rea- 
sons, only the R > 200 arcsec portion of the high-resolution 
Hi velocity field was used in the harmonic decomposition. 
The portion of the H i disc that we analyse using RESWRI is 
delimited by the dashed ellipses shown in Fig. [1] 



4.2.2 Results 

Figure [6] shows the results of the harmonic decomposition. 
The top three panels show the position angles, inclinations 
and systemic velocities from the circular tilted ring model 
fitted to the velocity field by RESWRI. These fitted param- 
eters agree well with those of the tilted ring model fitted 
to the velocity field by E2010a. They are also consistent 
with the corresponding optimal parameters determined by 
VELFIT (as listed in Table [1]). RESWRI creates a second resid- 
ual field by subtracting the fitted harmonic components of 
each ring from the observed velocity field. This residual field 
therefore provides a measure of the signal that was not cap- 
tured by the harmonic expansion. For all galactocentric radii 
larger than 200 arcsec we find a mean absolute residual of 
3.6 km s~^, suggesting that an harmonic expansion up to 
third ord er is capable of ca pturing most of the non-circular 
motions. ITrachternach et al.l (j2008|) find an average value 
of 2.9 km for the median absolute residuals of the 19 
THINGS galaxies in their sample. The average amplitude of 
the non-circular velocity components for their entire sample 
is 6.7±5.9 km s~^. The corresponding average for NGC 2915 
is 11.7 ± 4.0 km s"\ 

Within the theoretic a l fram ework of lFranx et al.1 (|l994h 
and ISchoenmakers et"aL Co and ci correspond re- 

spectively to the systemic and circular velocities of each ring, 
while the higher order terms measure the non-circular ve- 
locity components. A cursory examination of the harmonic 
decomposition results shows that the cs term is approxi- 
mately z ero over the full radi al range indicating, according 
to ISchoenmakers et al.l (Il997), that the corr e ct inc linations 
were fitted by RESWRI. ISchoenmakers et"al ] (|l997f ) do cau- 
tion, however, that in the case of an m = 2 spiral arm, 
this inclination is not necessarily equal to the true inclina- 
tion. However, in the case of NGC 2915, all of our three in- 
dependent disc inclination estimates, including the RESWRI 
estimate, are consistent with one another. Also shown in 
Fig. [6jl is the corresponding Vt profile of the optimal radial 
flow model fitted by velfit to the low-resolution H i velocity 
field. 



4.2.3 Discussion 

Perhaps the most striking result from the harmonic decom- 
position is the large amount of power seen in the si com- 
ponent (si > 10 km s"^ for R < 400 arcsec, and then 
5 km s~^ < si < 10 km s~^ for R > 400 arcsec). A strong 
si harmonic in the velocity field suggests the presence of an 
intrinsic m — 2 velocity distortion, often associated with a 
tri- axial halo. The problem, however, is that m — 2 spiral 
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Figure 6. RESWRI harmonic decomposition results (filled circles). Rows 1 and 2 show the circular tilted ring model parameters fitted 
by RESWRI to the velocity field. Rows 3-5 show the amplitudes of the harmonic components fitted to the residual velocity. Error bars 
represent the formal least-squares errors. In each panel of rows 4 and 5, the dashed line marks km s"-*^. The harmonic decomposition 
was only carried out for radii R ^ 200 arcsec. For comparison, the corresponding Vt profile of the optimal radial fiow model fitted by 
VELFIT to the low-resolution H I velocity field is shown in panel (d) as open circles. 



structure, which is clearly present in the Hi distribution, 
will also res ult in a sign i ficant si term in the observed ve- 
locity fi eld. I Wong et alJ (l2004l\. using the linearised equa- 
tions of Canzian AllenI ITqQtII for the velocity perturba- 
tions caused by a two- armed spiral density wave, point out 
that a characteristic kinematic signature in the velocity field 
due to such a perturbation is a sinusoidal variation with ra- 
dius of the si and S3 harmonic components. While it may 
be argued that some sinusoidal variation is present in the 
si term, the S3 term is close to zero for all radii, thereby 
providing evidence against the perturbational effects of a 
two- armed spiral density wave to the observed velocity field. 
Furthermore, the fact that the bi-symmetric models fitted 
by VELFIT to the velocity field are insensitive to spiral struc- 



ture, yet still suggest the presence of sizeable non-circular 
V2,t and V2,r velocity components, implies that the strong si 
term found by RESWRI cannot solely be due to spiral struc- 
ture. 

The strong si term could also be linked to an ax- 
isymmetric radial flow. In their description of the kine- 
matic signatures o f different types of non-circular motions, 
IWong et al state that a signiflcant si term together 

with a negligible S3 term can be caused by radial gas in- 
flow. We in deed see this behav iour in the harmonic terms of 
NGC 2915. IWong et al.l (|20o3 ) go on to list \ds3/dsi\ < 0.1 
over a significant radial range as a further kinematic signa- 
ture of radial flow. We measure an average absolute gradient 
of \dsz/dsi\ ^ 5.8 X 10"^ for the outer disc of NGC 2915, 
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consistent with the radial flow scenario. A further kinematic 
signature of radial flow is non-orthogonality of the kinematic 
major and minor axes. From the circular tilted ring model 
fltted by RESWRI to the H i velocity fleld, the position angle 
of the kinematic major axis is estimated to be ^ 287°. The 
iso- velocity contours in the H i velocity field suggest the posi- 
tion angle of the kinematic minor axis to be ^ 212° . The two 
kinematic axes thus appear to be non-orthogonal, differing 
in position angle by only ^ 75° , thereby further supporting 
the possibility of radial fiow. This evidence for radial fiows is 
further supported by the fact that velfit finds an axisym- 
metric radial fiow model to better match the observed Hi 
velocity field of NGC 2915 than does a bi-symmetric fiow 
model. This is true regardless of whether the bar angle in 
the bi-symmetric fiow model is kept fixed or is allowed to 
freely vary. 

IWong et al ] (|2004l l mention that if the si term is in- 
terpreted as evidence for radial fiow, then in a counter- 
clockwise rotating galaxy, a positive si term implies outfiow. 
Assuming that a galaxy's spiral arms are always trailing, 
NGC 2915 is rotating counter-clockwise, thereby implying 
the strong si term to be indicative of a radial outfiow of 
th e order of 5 — 17 km s~^. This result is similar to that 
of lGentile etHI (|2007l l who inferred the presence of radial 
infiows of the order of 5 — 13 km s~^ in t he Hi disc of 
the nearby dwarf irregular galaxy NGC 3741. Sancisi et ajj 
(,2008 ) hypothesise the existence of NGC 2915's extended H I 
disc to be as a result of the accumulation of in-fallen pris- 
tine gas from the surrounding inter-galactic medium. This 
gas, they say, forms "a reservoir of fresh gas for fuelling star 
formation in the inner regions" . The physical process of gas 
moving from the outer to the inner disc would, however, be 
associated with an axisymmetric radial inflow, and is hence 
difificult to reco ncile with our kinematic evidence for a ra- 
dial outfiow. We rk et al.l (2010) have recently reported on 
lower- and higher-than-expected oxygen abundances for the 
respective inner and outer portions of NGC 2915's Hi disc. 
The authors showed that the few isolated regions of low- level 
star formation within the outer disc cannot account for the 
measured oxygen abundance. They propose a metal-mixing 
scenario in which metals produced by the star forming core 
are radially redistributed towards the outer disc. Such a sce- 
nario is consistent with our kinematic evide nce for radial 
out fiows. F urthermore, several authors (e.g., 'Burea u et al.l 
[1999 : Mas set Bureaul [2003) strongly advocate the pres- 
ence of a bar within NGC 2915's co-rotation radius. If mass 
is able to fiow along this bar towards the centre of the galaxy, 
conservation of angular momentum dictates that some mass 
must also be sent outwards, away from the centre of the 
system. Our inferred outward radial motions are consistent 
with this physical scenario which could perhaps go some way 
in explaining the existence of NGC 2915's very extended H i 
disc. 

A major fiaw in the radial fiow argument is that the 
velocity field of a stationary bar potential re sembles ax- 
i symm etric infiow or outfiow. This, according to I Wong et al.l 
('2004'), is because for a fiat rotation curve the S3 and C3 
terms are expected to vanish in the case of a stationary 
pertu rbation to the potential (demonstrated bv [FYanx et al.l 
I1994I ). Bars are indeed thought to be the most effective 
mechanism for transporting gas to the centre of a galaxy 
dCombed 1 19991) , and so the kinematics alone cannot dis- 



Table 2. Summary of Hi velocity field harmonic decomposition 
results for NGC 2915. 



Perturber 
type 



Evidence 
for 



Evidence 
against 



m = 2 spiral 
density 



Strong si 
term 



Tri- axial 
halo 



Axisymmetric 
radial flow 



Strong si 
term 



Strong si term 
and weak S3 term 

\ds3/dsi\ ^ 

Non-orthogonal 
kinematic axes. 

Well-fltted VELFIT 
axisymmetric flow 
models. 



S3 term does 

not vary 
sinusoidally. 

Well-fltted VELFIT 
bi-symmetric flow models 
are insensitive to 
spiral structure. 



Contributions from 
radial flow, and 

m = 2 perturbation 
due to H I bar. 



Indistinguishable from 
elliptical streaming 
in a stationary bar 
potential. 



tinguish between t he stationary bar potential and the ra- 
dial fiow scenarios. I Wong et al.l (|2004l ) carried out harmonic 
decompositions of the CO and Hi velocity fields of seven 
nearby spiral galaxies in order to search for evidence of ra- 
dial gas fiows. They found no unambiguous evidence for ra- 
dial fiows, concluding that "the inherent non-axisymmetry 
of spiral galaxies is the greatest limitation to the direct de- 
tection of radial fiows". Thus, although the si term pro- 
vides very suggestive evidence for gas outfiow, it should not 
be treated as conclusive evidence for such gas kinematics. 
Elliptical streaming of the gas in a bar-like potential can 
still contribute significantly to the observed non-circular mo- 
tions. The various perturbations to the potential described 
in this section are summarised in Table [21 together with a 
list of evidence for and against each perturbation type. 



4-2.4 Elongation of the potential 

NGC 2915 is known to be extremely dark-matter- 
dominated, with a total mass to B-band light ratio of 
^ 140 M0/L0 (E2010a). If the non-circular components 
measured in the velocity field are treated as being due to 
a tri- axial dark matter halo, th ey can be used to estimate 
the axial ratio of the potential. ISchoenmakers et al.l (|l997l l 
present an expression for the global elongation of the poten- 
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tial which depends on the si and S3 harmonics measured by 
RESWRI. In the case of NGC 2915, however, the observable 
signatures of a possible radial outflow and an m = 2 spiral 
density wave perturbation to the potential will couple signif- 
icantly with those of an elongated potential, the axial ratio 
for which therefore cannot be confidently qua ntified using 



tor wnicn tnererore cannot be connaent i y qua n 
the expression from ISchoenmakers et al. (Il997l l. 

As exp lained i n ISpekkens SellwoodI |2007) and 
ISellwood Zanmar Sanched (I2OIOI). the components re- 
turned from VELFIT bi-symmetric flow models are not 
sensitive to these confounding factors. For a mildly dis- 
torted, slowly rotating potential and a flat rotation curve, 
ISellwood Zanmar Sanchez! (|2Q1Q ) estimate the axial ratio 
of the potential as 



Vt 



Vt + 2V2, 



1/2 



(7) 



We use the Vt and V2,r profiles from our bi-symmetric flow 
models fitted to the low-resolution Hi velocity field to es- 
timate qcf). We do this for the model in which = 322° 
was kept fixed, and for the model in which it was allowed to 
freely vary (attaining a best-fitting value of ^J, = 319.2°). 

The two models yield very similar radial profiles, 
shown in Fig. [71 In both cases, the axial ratio for the po- 
tential is close to unity for nearly all radii, generally being 
> 0.9. Both bi-symmetric flow models have slowly varying 
V2,t and V2,r profiles for 250 arcsec < < 400 arcsec. Addi- 
tionally, the amplitudes of the non-circular velocity compo- 
nents are less than 10 percent of the corresponding circular 
velocities. We therefore treat the estimates within this ra- 
dial range as being the most reliable. Average axial ratios of 
qcf) = 0.92zb0.21 and = 0.92zb0.35 are obtained within this 
radial range for the models in which ^J, is fixed and allowed 
to vary respectively. These estimates are consistent with an 
axial ratio of = 1, and are consistent with the estimated 
lower limit of > 0.98 from ISellwood Zanmar Sanched 
(|20ld l for the axial ratio of NGC 3198. We stress, however, 
that Eqn.[71m// not yield reliable estimates when the po- 
tential is significantly distorted or when it is fast- rotating. 
Neither of these possibilities can be confidently ruled out 
for the case of NGC 2915. While our results are strongly 
suggestive of a perfectly round halo for NGC 2915, other 
possibilities cannot be excluded due to the limitations of 
the method discussed above. 



5 SUMMARY AND CONCLUSIONS 

We have searched for non-circular flows within the extended 
Hi disc of NGC 2915 using two independent methods. The 
resulting non-circular velocity components have been inter- 
preted in the context of axisymmetric and non-axisymmetric 
perturbations to the system's gravitational potential. The 
results allow us to probe the gravitational potential of the 
system out to radii far beyond the edge of the stellar disc, 
into the regions of the galaxy that are completely dark- 
matter-dominated. 

The first implemented routine, velfit, fits a specified 
kinematic model to the entire Hi velocity field. Models were 
fitted to both high- and low-resolution versions of the obser- 
vations, the corresponding results of which are consistent 
with one another. Axisymmetric radial flow models best 
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Figure 7. Radial distribution of the axial ratio of the potential 
of NGC 2915. Closed and open circles correspond to our fitted 
bi-symetric flow models in which cf)'^ was held fixed and allowed 
to freely vary respectively. For the sake of clarity, the open circles 
have been shifted by —5 arcsec. Error bars were estimated using 
Gaussian error propagation of the non-parametric estimates of 
the true uncertainties in Vt and V2,r, generated by applying a 
bootstrap method. 



match the observed Hi velocity field. The optimal radial 
flow models suggest the presence of radial flows as large 
as 20 km s~^ over a large fraction of the Hi disc. We ar- 
gue that, at least at inner radii, these inferred radial flows 
could be the result of the deposition of kinetic energy into 
the inter-stellar medium by massive stars and supernovae 
explosions, velfit was also used to fit bi-symmetric flow 
models to the H I velocity field. The optimal models suggest 
bi-symmetric flows of ^ 8-20 km s~^ throughout most of 
the H I disc. These flows occur along an axis oriented ~ 18°- 
30° relative to the major axis of the disc, thereby linking 
them to the central bar-like feature seen in the mass distri- 
bution of NGC 2915. These flows could therefore be caused 
by streaming gas motions within a bar-like perturbation of 
the potential. 

NGC 2915 is extremely dark-matter-dominated. The Vt 
and V2,r profiles from our fitted bi-symmetric flow models 
have been used to estimate the axial ratio of the potential 
under the assumption that it is mildly distorted and slow- 
rotating. Our results are consistent with the potential being 
perfectly round, with an average axial ratio of ^ 0.92. We 
stress, however, that because of the underlying assumptions 
used to obtain this result, it cannot be treated as indis- 
putable evidence against a possibly aspherical halo. 

The second routine implemented to quantify the non- 
circular motions is RESWRI. Evidence is found for an m = 2 
term in the Hi velocity field. The possibility of elliptical 
streaming in a stationary bar-like potential detracts from 
the otherwise good evidence that is provided by the strong 
and weak si and S3 terms for an axisymmetric radial outflow 
of the order of 5-17 km s~^. Such outflows, ^ 6-20 percent 
of the outer Hi disc rotation velocity, do not provide intu- 
itive insight into the formation history of the Hi disc if it 
is thought of as having accumulated from in-fallen pristine 
gas from the surrounding inter-galactic medium. Under such 
circumstances, the in-fallen gas would be associated with a 
radial inflow rather than an outflow as it is transported to 
the central regions of the galaxy. It is possible, however, that 
some material is being re-distributed into the outer disc in 
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order to conserve angular momentum as other material flows 
inwards along a bar. A clear understanding of the formation 
history of NGC 2915's extended Hi disc is nonetheless still 
elusive. The limitation on the search for indirect evidence 
of radial flows set by the inherent non-axisymmetry of the 
galaxy gives rise to the need for deeper Hi observations of 
NGC 2915's outer disc in order to directly detect possible 
cold gas accretion processes and to carefully search its local 
environment. 
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